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PREFACE 


» 


This report describes part of a continuing research program in remote 
sensing of the environment. The research is being conducted for NASA’s 
Lyndon B, Johnson Space Center, Houston, Texas, by the Environmental 
Research Institute of Michigan (ERIM) , formerly the Willow Sun Laboratories, 
Institute of Science and Technology, The University of Michigan. The main 
objective of this program is to develop remote sensing as a practical tool 
for obtaining extensive environmental information quickly and economically. 

Remote sensing of the environment involves the transfer of radiation 
from a target on the Earth’s surface through the atmosphere to a sensor 
located at a point either within or above the atmosphere. The atmospheric 
medium can alter the radiation by absorption and scattering. In this report 
we examine the effects due to the spatial variation of the atmosphere and 
also due to influence of background on target as a result of scattering. 
Results arc presented in terms of classification accuracy for variable 
atmospheres and in terms of percent change in target radiance for the 
target-background interaction problem. 

The research described in this report was performed under Contract 
NAS 9—1412 3, Task VIII, and covers the period from 15 May 1974 through 
14 March 1975. Dr, Andrew Potter served as Technical Monitor. The program 
was directed by R. R. Legault, Vice President of ERIM, J. D. Erickson, Head 
of ERIM’s Information Systems, and Analysis Department, and R. F. Nalepka, 
Principal Investigator and Head of ERIM’s Multispectral Analysis Section. 

The ERIM number of this report is 109600-15-F. 
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and L. R. Ziegler. The author also thanks D. Dickerson and L. Parker for 
their secretarial assistance in the preparation of this report. 
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1 

SUMMARY 


There has been a continuing question in multispectral remote sensing 
data analysis as to the extent that the atmosphere affects the classifi- 
cation of various materials on the Earth's surface. In previous studies 
we have developed a radiative- transfer model and have used that model to 
examine the variation in simulated multispectral data. The model was a 
simplified yet practical model which included the effects of absorptive 
and non-absorptive haze conditions. The main limitation of the model, as 
far as remote sensing data analysis is concerned, was its one-dimensional 
nature, i.e., no spatial variations in the horizontal plane were considered. 

In the present treatment we have extended the model to include non- 
uniform atmospheres and have examined the influence of the resulting changes 
in radiance on the accuracy of classification of various classes of objects. 
Various atmospheric gradients of optical thickness were considered throughout 
a LAHDSAT frame and fields were produced artificially by random sampling 
from real data. Classification accuracy was then determined for each field 
and analyzed as a function of visual range. The results indicate that as 
the visual range atmospheric gradient increases the classification accuracy 
decreases. An atmospheric correction algorithm is also presented which can 
be used to eliminate the systematic variations attributable to the atmosphere. 

Potentially, a more important effect is the influence of background 
materials on the target as a result of atmospheric scattering. In particular, 
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we would like to know what effect the proximity of various background elements 
has on the spectral character of the target. If there is a strong influence, 
then the probability of misclassification could increase significantly. A 
detailed three-dimensional treatment of this problem has been carried out 
and the influence of background on target has been analysed for various 
environmental conditions. In particular, we found that for dark targets 
and bright backgrounds the per cent change in the target radiance can be 
quite large, i.e,, 50 “ 70 %. We also examined this effect in terms of sun 
angle and visibility and found that the maximum variability occurs in the 
normal range of operat tonal parameters . The gradual change in target radiance 
is given in terms of varying amounts of background material, thus simulating 
a variety of possible geometric background field configurations . An algorithm 
which can he used to remove this effect from actual multispectral data is also 
described . 


2 

INTRODUCTION 

In remote sensing of features on Earth's surface the min purpose 
of multispectral data analysis is to discriminate among various materials 
j 3 nrl to recognize classes of objects. There are many effects which must 
be included in such investigations . There is never a 100 per cent c las si- 
fication because of various noise factors in the system. Some of these 
factors are: intrinsic variations within a class at a particular time, 

such as the natural variability in a wheat field at a particular point in 
the growing season; mixtures of various materials in the same area; intrinsic 
temporal variations of a material; intrinsic goniometric properties; cli- 
matic and weather changes; atmospheric effects; and detector noise. Many 
of these effects can be dealt with in a statistical manner, and using 
statistical decision theory one can arrive at a signature or class identi- 
fication. 

The remote sensing of terrain features means that we are necessarily 
receiving radiation which has passed through part or all of the Earth's 
atmosphere. As a result, it has long been suspected that the atmosphere 
c,nn cause significant changes in signatures and thereby reduce the accuracy 
of classification. In the current investigation, we shall consider primarily 
two atmospheric effects; changes in classification accuracy for spatially 
varying atmospheres, and the target-background interaction, i.e., the 
influence of background elements on the target as a result of atmospheric 
scattering . 
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A rigorous mathematical treatment of radiative-transfer in a spatially 
varying atmosphere is an intriguing hut extremely difficult problem. 
Nevertheless, under certain circumstances, simplifying assumptions can 
be made and a reasonably accurate practical solution can be found. In 
this study we shall analyze various atmospheric conditions using a simpli- 
fied radiative-transfer model to simulate multispectral data for different 
classes of materials. We can then vary the atmospheric state and see what 
the resulting change is in the classification accuracy. 

It turns out that the second atmospheric effect, i.e., the target- 
background interaction is in many cases more important as far as changes 
in spectral radiances are concerned, than variable atmospheres. We want 
to know for example, the change in spectral radiance for a target surrounded 
by various patterns of background materials . Also, we want to know under 
what conditions this effect is most important. 

In the mathematical derivations of the basic equations very few approxi- 
mations are made. The most important ones are: 1) a uniform atmosphere, and 

2 ) single scattering. Multiple scattering is included in the determination 
of "pure’* atmospheric path radiance and in the downward diffuse irradiance. 

A somewhat rigorous treatment is necessary because we do not want to intro- 
duce anomalous mathematical uncertainties . We do want whatever variations 
occur to reflect intrinsic changes in the environmental conditions and not 
be the result of mathematical approximations . In the actual analysis we have 
made other assumptions however, in order to simplify the computational analysis. 
For example, we assume a checkerboard pattern of uniform square fields each 
being a Lambert reflector. The general mathematical model, however, is capable 
of using any complex pattern and a surface with any goniometrie property. 
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A further question arises; can we develop an algorithm so as to remove 
this effect from real data? The answer is yes s providing that we have 
some knowledge of the atmosphere and the surface pattern of fields. In 
the last section we devise an algorithm which can be used in the analysis 
of actual multispectral data. 

Having solved these problems it would now be appropriate to consider 
both the variable atmospheric state and the interaction problem simultaneously 
in order to determine classification accuracy. 
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VARIABLE ATMOSPHERES 

In this section we shall examine the variation of radiometric quantities 
in terms of the parameters used in the radiative-transfer model. Making 
use of these parameters we will then construct hypothetical variable atmospheric 
states for which we will implement the atmospheric correction algorithm. The 
percent classification will be calculated for the various atmospheric conditions. 
3.1 OPTICAL PARAMETERS 

In order to construct variable atmospheres we will first define the basic 
optical parameters which are used in the radiative- transfer model. 

3.1.1 ATTENUATION COEFFICIENTS 

The attenuation of a collimated beam of radiation is expressed by the 
following equation? 

I^(x) » X^(o) e * (1) 

when I-^(x) and (o) are the spectral intensities of the beam at distance x and 
at the origin respectively and is called the spectral extinction coefficient. 
This is a measure of the total removal of radiation from the initial beam as 
a result of both scattering and absorption in the medium. 

In 1924, Koschmieder [1] defined a distance (visual range's such that the 
limiting contrast between an object and its background (pure atmosphere) is 0.02. 
This definition provides us with the following useful formula which relates 


Koschmieder, Beitr, Phys. freien Atm., Vol. 12, 19,24 pp. 33-53, 

171-181. 
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visual range to the volume extinction coefficient at a wavelength of 0.55 pm: 


V(km) = (2) 

(km x ) 
o 

where X = 0.55, This relationship holds ever a wide range of distances 
and is illustrated in Fig, 1. 

For realistic atmospheres we must consider not only the gaseous com- 
ponent but also the semi-permanent suspension of particulates t referred to 
as the aerosol component. Thus, we have 

K^(s) - K^(s) + K^(2> ^ 

R A 

where k^(z) is the well-known Rayleigh coefficient at altitude s and tc^(z) 

is the corresponding aerosol coefficient. Besides the volume extinction 

coefficients we also have volume absorption and volume scattering 

coef f icient s , i . e . , 

k^(z) = “ x (z) + B x (s) (4) 


where the absorption coefficient a^(z) aad the scattering coefficient 3^(z) 

for the gaseous component is known by experiment and the corresponding 

aerosol coefficients are determined by Mie theory, i.e., 

00 

a^(z) = | n(r,s) o^(r)dr (5) 

o 


and 




CD 

n(r s z) cr!f(r)dr 
A 

o 


<6) 
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where n(r,s) is the particle size distribution function for particles of 
radius r and o^(r) and o^(r) arc the absorption and scattering cross 
sections for uniform spherical particles. Details of the properties of 
these coefficients were determined for realistic atmospheres by Turner [Z] . 


3.1.2 OPTICAL DEPTH 

A more convenient parameter to use in radiative- transfer theory than 
the attenuation coefficient is the optical depth t. This is a dimensionless 
quantity which is a measure of the turbidity of the atmosphere. Mathematically, 


it is expressed as follows: 



rCO 

K, (z)dz 


(7) 


where is the spectral optical depth at altitude h. Using this definition, 
the optical depth of the entire atmosphere, called the optical thickness t q ^, 
is given by 

t oA - j K A <z)dz (8) 

0 

Typical values of in the visible part of the electromagnetic spectrum 

range from ^ 0.03 for clear atmospheres to ^ 3 for very hazy atmospheres. 

From a practical standpoint, the spectral optical thickness is an 

easily measurable quantity. Using a spectrometer or soma multispectral 

radiometric device one can measure the spectral ir radiance of the direct 

-1 

sunlight, EjjCA) for a specific solar zenith angle 0 Q (=cos ji Q ) and determine 


F21 

1 R, E. Turner, Radiative Transxer in Seal Atmospheres, Report Ho. 
190100-24— T, Environmental Research Institute of Michigan, Ann Arbor, 1974. 
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the optical thickness from the following formulas 


P E (X) 

n O O 

T oX " P o E,(X) 
a 


( 9 ) 


where E q (X) is the known extraterrestrial value of the solar irradiance. 
Making use of data obtained from light scattering experiments in the 
atmosphere, Elterman [3J constructed model atmospheres for various hazes. 
His data relate the optical depth for any altitude to the wavelength and 
visibility. 

3.1.3 SCATTERING PHASE FUNCTIONS 


Another quantity of importance in radiative- transfer analysis is the 
single-scattering phase function £)^(r s 0,$) . It describes the angular 
scattering characteristics of a particle of radius r and can be thought of 
as a probability density. It is normalised such that the integral of the 
phase' function over all angles is unity, i.e.. 


I I d ^ cos6 ^ d ^ “ 1 ( 10 > 

For particles which are small as compared to the wavelength of the incident 

radiation the angular variation in f3^{r s 0,4>) is almost like that for dipole 

o 

scattering, i.e., ^ 1 + cos X» where X is the scattering angle. For large 
particles, however, the scattering is predominantly in the forward direction. 
Deirmendj ian [^ ] has calculated the phase function integrated over various 
particle size distribution in order to obtain the so-called polydisperse 

_ 

1 L. Elterman, Vertical -Attenuation Model with Eight Surface Meteorological 
Ranges 2 to 13 Kilometers, Report No. AFCRL-70-02GQ, Air Force Cambridge Research 
Laboratories, Office of Aerospace Research, Bedford, Mass., 1970. 

[4] 

D . Deirmendjian, Electromagnetic Scattering on Spherical Polydispersions, 
Elsevier, New York, 1969, p. 78. 
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scattering phase function, V--. have also run computer programs to determine 
the scattering phase functions and attenuation coefficients for various 
atmospheric states characterized by particle sizes and compositions, A 
detailed treatment of these optical parameters is beyond the scope of this 
report, but can be found in the author's previous report [ 2 ], 

3.2 RADIOMETRIC QUANTITIES 

The general solution of the radiative-transfer equation is the following: 


L = LT + L 
o p 


< 11 ) 


where L q is the spectral radiance at the surface, T is the spectral transmittance 
from the surface to the sensor, and is the spectral path radiance, i.e., 
that radiance which results from all multiply-scattered radiation along 
the path from the target to the sensor. The radiative- transfer model developed 
at ERIM allows one to calculate all of the quantities in Eq. (11) for various 

values of the independent variables on which the radiometric quantities 
depend. The explicit dependence is as follows: 

L - L(A,T o , z, 0, tf>, 6 q , p, p) (12) 


1* 0 k Q (A, Tq, 0 , <}>, 0 q, p, p) 


(13) 


T = T(A, t q , z, 8 ) 


(14) 


l*p ” ^(^9 Iq 9 z ® G» <f>> P * p) 


(15) 


m 

1 R. E. Turner, Radiative Transfer in Real Atmospheres, Report No. 
1901G0-24-T, Environmental Research Institute of Michigan, Ann Arbor, 1974. 
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where X is the wavelength, t q the optical thickness, 6 the nadir view angle, 

$ the azimuthal angle between the view plane and solar plane, 8 the solar 

o 

zenith angle, z the altitude of the sensor, p the target reflectance, and 
p the background albedo. In addition, for the case of perfectly diffuse 
(Lambertian) surfaces the surface radiance is given by 



where 

E = E(X, t o , 0 O , P, P) (17) 

is the total downward spectral irradiance on the target =, By specifying 
various input parameters the radiometric quantities can be calculated. 

Figure 2 illustrates the input-output characteristics of the basic radiative- 
transfer model developed at ERIM. 

Of all radiometric quantities the easiest to measure is probably the 
direct solar irradiance. This measurement gives us the optical thickness of 
the atmosphere, which in turn allows us to calculate the transmittance, I.e., 


T 

o 



where T is the vertical transmittance of the entire atmosphere. The path 
o 

radiance is vsually more difficult to measure accurately. The ideal way 
to determine it is to measure the radiance of an absolutely black infinitesimal 
target. Then, since the reflectance p is zero, we have, according to Eq. (11) 

L » L p (18) 

Although there do exist surface materials which have very low reflectances, 
there is still some significant amount of radiation from the target. 

In order to understand the relationship of these various quantities we 
shall plot a series of graphs. 
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Solar Zenith Angle (Time,Date) 
Zenith View Angle 
Nadir View Angle 
Relative Azimuthal Angle 
Target Reflectance 
Background Reflectance 
Target Radius 


Bate 

Wavelength 

Base Atmospheric Temperature 
Base Atmospheric Treasure 
Visual Range (Optical Thickness) 
Aerosol Refractive Index 
Aerosol Size Distribution 
Sensor Altitude 



Air Refractive Index 
Single-Scattering Albedo 
Single-Scattering Phase Function 
Optical Depth 


Transmittance 
Ir radiance 
Surface Radiance 
Beam Radiance 

Singly-Scattered Solar Radiance 
Doubly-Scattered Solar Radiance 
Singly-Scattered Surface Radiance 
Path Radiance (Multiple Scattering) 
Sky Radiance (Multiple Scattering) 
Total (Beam plus Path) Radiance 


FIGURE 2. ERB1 RADIATIVE TRANSFER MODEL 
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figure 3 illustrates the radiation of the direct solar and downward 

diffuse irradiances as a function of the atmospheric optical thickness and 

surface albedo. The direct irradiance decreases exponentially and is 

independent of surface conditions since it just represents the amount of 

radiation, which is not lost from the solar input. For an optical thickness 

of zero (no atmosphere) there is, by definition, no diffuse component, hut 

as x increases, some of the energy lost by the direct beam goes into the 

o 

diffuse downward component. It is interesting to note, however, that a white 
surface only raises the diffuse downward irradiance by about 30% over that for 
a black surface in the case of a thick atmosphere. For a thin atmosphere, 
however, the white surface produces an irradiance almost three times that for 
a black surface. This effect is illustrated in Fig. 4. The physical 
principles can be understood better by looking at the total downward irradiance 
as a function of transmittance T . This variation is depicted in Fig. 5, For 
darker surfaces the irradiance decreases with increasing turbidity because 
less radiation is able to reach the surface. On the other hand, for brighter 
surfaces, more radiation is reflected and is essentially trapped under the 
ever thickening turbid atmosphere. 

All of the radiometric quantities are interrelated, sometimes in a 
rather complicated way. Figure 6 illustrates the dependence of both 
radiance on irradiance for various transmit tances and surface albedos. 

It is interesting to note that both of these quantities depend upon surface 
albedo linearly. Also, for larger values of transmittance we see a very 
large change in path radiance but a corresponding very small change in 
irradiance . 
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FIGURE 3 


Dependence of Direct Solar and 
Diffuse Downward Irradiance on 
Optical Thickness and Surface 
Albedo. Wavelength ■ 0.55 ym, 
Solar Zenith Angle * O'*. 
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FIGURE 4 


snce of Diffuse Downward 
on Surface Albedo for 
.lv thin and thick Atmospheres 
»gth - 0.55 ym, Solar Zenith 
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FIGURE 5 

TOTAL DOWNWARD IRRAD LANCE AS A 
FUNCTION OF TRANSMITTANCE FOR 
VARIOUS SURFACES. 

SOLAR ZENITH ANGLE - 30°, 
WAVELENGTH - 0.55 ym 
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FIGURE 6 


DEPENDENCE OF PATH RADIANCE ON TOTAL 
IRRADIANCE FOR VARIOUS TRANSMITTANCES 
AND SURFACE ALBEDOS. SOLAR ZENITH 
ANGLE = 30°, NADIR VIEW ANGLE - 0°, 
WAVELENGTH = 0.55 pm 
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For practical cases of fcransmittances which are usually encountered 
we can relate them to visual ranges. This is shown in Fig. 7. 

We should now examine the total radiance. If we let the background 
albedo be held constant and vary the target reflectance, then the total 
radiance increases in a linear way as depicted in Fig. 8. Likewise, holding 
the target reflectance constant, we can see the effect of a variation in the 
background albedo. This effect is illustrated in Fig. 9 for four atmospheric 
states. It should be noted from both of these graphs that there exists a 
wide range of atmospheric states and a very small range of reflectances which 
can give approximately the same total radiance, i.e., given a certain 
reflectance, we can have almost any atmospheric state and the total radiance 
remains the same. 

We can ex amin e the dependence of total radiance on the irradiance. This 
relationship is illustrated in Fig. 10. The dependence of each on reflectance 
is linear but for specific transmittances the dependence is definitely non- 
linear. For example, for a surface with a reflectance of 50% the radiance 
decreases with increasing transmittance up to a point, but then the radiance 
increases rapidly with a small increase in transmittance. If, however, the 
surface reflectance is very high, say 100%, then the radiance increases as the 
transmittance increases. These effects are due to the fact that darker surfaces 
will get darker as the haze decreases but brighter surfaces will get even 
brighter. It should also be noted that if the total radiance is measured 
from a sensor aboard a satellite and the transmittance is determined, then a 
unique value for the reflectance exists. This is true however, only for p= p. 
Nevertheless, the same procedure can be followed for more general conditions 
in which p ^ p. It is a mathematically well-determined problem since the 
irradiance is also a measurable quantity. 
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FIGURE 8 


Dependence of Total Spectral 
Radiance on Target Reflectance 
and Optical Thickness for a 
Background Albedo of 0.4. 
Wavelength “ 0.55 urn, 

Solar Zenith Angle “ 60° 

Nadir Angle >0° 


TARGET REFLECTANCE p 












FIGURE 9 


Dependence of Total Spectral 
Radiance on Surface Albedo and 
Optical Thickness for a Target 
Reflectance of 0,2, 

Wavelength - 0.55 ym. 

Solar Zenith Angle » 60° 

Nadir Angle » 0° 
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FIGURE 10 

DEPENDENCE OF TOTAL RADIANCE ON TOTAL 
IRRAD LANCE FOR VARIOUS TRANSMITTANCES 
AND SURFACE ALBEDOS. SOLAR ZENITH 
ANGLE = 30°, NADIR VIEW ANGLE = 0°, 
WAVELENGTH = 0.55 pm 
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Finally, we can analyze the relationship between total radiance and 
path radiance (Fig. 11). Thus, for dark surfaces a hazy atmosphere produces 
a higher total radiance and path radiance than for clear atmospheres. For 
bright surfaces, however, the total radiance for hazy atmospheres is less 
than that for clear atmospheres. This effect is demonstrated by looking 
at ERTS photography (Turner et al., [6 ]). 

It is assumed that in any analysis of multispectral data the wavelength, 

altitude, and geometric conditions are known. The unknown variables are 

optical thickness t , target reflectance p, and background albedo p. The 

o 

values of these variables can be determined, however, by simultaneous 
measurement of the quantities L, L^, T, and E. Thus, for uniform surfaces 
the state of the atmosphere and surface reflectances can be determined by 
analyzing multispectral data and auxiliary data. It is not always practical 
to measure the auxiliary data, however, and also the usual surface pattern is 
not uniform. The more complex case of a non-uniform surface will be treated 
in a later section. 

3.3 CORRECTION ALGORITHM 

We shall now consider a method which will allow us to eliminate the 
variable effects of the atmosphere. For simplicity, we shall investigate 
the case of satellite remote sensing, i.e., at an altitude where the optical 
depth is zero. The basic equation for total radiance is then 

L(o, y, <J>) = L (t , y, 4>)T(o,y) + L (o, y , 4>) (19) 

o o p 

^R. E. Turner, W. A. Malila, R. F. Nalepka, F. J. Thomson, Influence 
of the Atmosphere on Remotely Sensed Data, Proceedings of Society of Photo- 
Optical Instrumentation Engineers, 1975. 
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FIGURE 11 

DEPENDENCE OF TOTAL RADIANCE ON PATH 
RADIANCE FOR VARIOUS TRANSMITTANCES 
AND SURFACE ALBEDOS. SOLAR ZENITH 
ANGLE = 30°, NADIR VIEW ANGLE = 0°, 
WAVELENGTH = 0.55 urn 
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where the zero means x = o. The surface radiance is given by 


2 7T 2 


W M) 


y , p(y,*,u’,*’)tL sun (t o ,-n , ,t’) + L (T o ,-^ , ,* , ))d u , d*' 


o o 


where the direct solar radiance is 

-x /y 


L (x ,-y,^) = E e S(y-y ) ) 

sun o o o o 

and the second part of Eq. (20) is defined as the diffuse sky radiance 

Hence, we have 

-t /y 

L (x , y , 4>) * y E e ° ° p(y,<j>,-y ,* ) 

o O 0 0 o o 


( 21 ) 


( 22 ) 


for a general surface. We can now define an intrinsic radiance L^.(y J <)>) 
as being that radiance which we would get if no atmosphere were present, i.e. 


Lj (y,(j>) = y Q E o p(y,<}>,-y o ,<|> o ) 


(23) 


Thus, Eq. (19) becomes 

L(o,y s $) = L I (y,<^)T(o,y o )T(o,y) + L n (y,<t>)T(o,y) + L p (o,y,<|>) (24) 

The intrinsic radiance can therefore be put in the form 

L_(y,<f>) = R (y,y ) L(o,y,4>) - S (y,y ,<|>) (25) 

I o o 

where Lj(y,d>) is the "corrected" radiance and L(o»y,ij>) is the measured value 
obtained by remote sensing. The correction functions R(y,y^) and S(y,y o ,4>) 
are either determined experimentally or by model calculations. If we assume 
Lambertian surfaces then the correction functions can be calculated easily. 
In this case 

L d (y,<|>) = ^ E_ (26) 


( 20 ) 
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where E_ is the downward diffuse irradiance on the target. The intrinsic 
radiance in this case becomes 

Lt(p»<J>) « -r p E ^ 27 ) 

■L TT O O 

Now, we also have 

E - E - p E T(p ) (28) 

- - o o o 

where E_ is the total downward irradiance on the target. Using these 
relationships on Eq. (24) we get 

R(y,u ) = ; 

T(u o )[T(y)-l3 + E_/y o E o (29) 

S(p,p 0 ,<|>) ■ R(y,n o ) L p (o,u,4>) (30) 

where all quantities in Eqs. (29) and (30) are easily measurable except 
the path radiance. 

Standard procedures were developed (Turner & Spencer [6], Horvath et.al., 

[7]) for the correction of raultispectral satellite data for atmospheric variations 
involving sun angle, surface albedo, and atmospheric tra^ — ittance. The methods 
are not as general as the procedure described above but do provide the user with 
a simple algorithm which can be used to quantify the atmospheric variations. 

3.4 VARIATION OF PARAMETERS 

Ir* order to understand the variation in multispectral data as a result 

of a variable atmosphere we will construct an atmospheric state and calculate 

radiance at various points of the scene. These will then represent the actual, 

^R.E. Turner & M.M. Spencer, Atmospheric Model for Correction of Spacecraft Data, 
Proceedings of the 8th Int. Conf. on Remote Sensing of Env., Ann Arbor, 1972 

^R. Horvath, M. Spencer, .& R. Turner, Atmospheric Correction & Simulation of Space 
Acquired Remote Sensor Data: 0. 4-1.0 ym Spectral Range, Report 10657-5-F, WRL, 
Institute of Science & Tech. Univ of Mich. 1972. 
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measured radiance for a specific atmospheric condition. Then* selecting a 
smaller number of points which would simulate auxiliary measurements of optical 
thickness, we use these data points to construct an approximate surface and 
then calculate new radiance. Comparing the actual with the calculated then 

provides us with a measure of the atmospheric variability. 

An we have seen, for a given visual range the optical thickness may 

vary around some mean value x . Let us therefore consider a 21 x 21 

o 

element frame of data with a spatially constant surface albedo and allow 
the optical thickness to have maxima and minima at N points in the scene, 
i.e., 

- J -[a i (x-x ± ) 2 + b,(y-y.) 2 ] 

? (x,y) - t + l a e 1 1 (31) 

i=l 

where the is the value of the increase or decrease in the mean optical 
thickness at the i^h point and the a^, b^ are the widths of the 

Gaussian distributions. Using the radiative-transfer model we can cal- 
culate the actual radiances for this surface at any point (x,y) . Then, 
choosing n sample points corresponding to optical thickness measurements 
we determine the corresponding T o (x^,yj), (j=l,..,n) and using regression 
analysis determine the best surface through the n points. As n ‘increases 
the radiances as calculated by using the regression surface should approach 
those determined by the exact calculation. 

Case I: Four Point Atmosphere 

As an illustration of this method let us consider the simplified 
case of circular contours at four points in a 21 x 21 frame. Thus 
3 i = ~ 3 f° r consider two visibilities with 

fluctuations in the optical thickness surface consistent with Eq. (33). 

The parameters are given in Table 1. 
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TABLE 1. PARAMETERS FOR MULTI-GAUSSIAN OPTICAL THICKNESS SURFACE 


i 


V = 10 

km (x - 0.611) 
o 


V = 23 km 

(t q = 0.379) 

i 

X. 

X 

y i 

°i 

X i 

y i 

°i 

1 

6 

6 

0.072 

6 

6 


2 

6 

16 

0.241 

6 

16 

0.125 

3 

16 

6 

-0.241 

16 

6 

-0.125 

4 

16 

16 

-0.072 

16 

16 

-0.038 


The smaller absolute quantities correspond to a 15 per cent fluctuation in 
the aerosol optical thickness and the larger absolute quantities refer to a 
50 per cent fluctuation. Using these parameters in Eq. (31.) we calculated 
an "actual" radiance L for each of the 441 points of the frame. Then 

cl 

selecting n(<441) points according to the arrays shown in Fig.. 12 we used 

regression analysis to determine the best surface through the n points. 

Using this surface we then calculated radiances L for each of the 441 points. 

c 



FIG. 12. OITICAL THICKNESS EXPERIMENTAL ARRAY 
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If the t q regression surface corresponded exactly to the known t q surface 

all the radiances should be the same. In reality there is a spread. For 

the 1, 9, and 25 sample cases we plot histograms in Fig. 13. As can be 

seen there is a considerable improvement as the number of "experimental" 

sample points increases. Figure 14 illustrates the difference between the 

two arrays involving five sensors. Figure 15a illustrates the spatial 

variation in the actual radiance based upon the actual t surface (Fig. 15b) 

o 

according to Eqn. 31. The large peak occurs in the lower left corner and 
the big valley is in the upper right corner. For n-1, i.e., only one 
sample point in the center of the frame the regression surface is a plane. 
Hence, a contour plot of L /L should merely reflect the variation of the 

C Si 

t q surface. This is demonstrated in Fig. 16a. For the two cases of five 
sensors, however, the ratio varies somewhat as is illustrated in Figs. 

17a and 17b. The regression surfaces in these cases are quadratic in 
the x direction but linear in y. For the nine sensor array there is 
some improvement as shown in Fig. 16b. Finally, the 25 sensor case is 
depicted in Fig. 18. Here the regression surface is a plane instead of 
a two-dimensional quadratic. This undoubtedly arises because of the 
symmetry introduced by our original choice of a multi-Gaussian surface. 
Hence, the improvement is less than in the nine sensor case. This means 
that one really should use regression analysis for much higher order 
algebraic surfaces or for multi-Gaussian regression surfaces. Such 
analysis is now continuing. As proof of the concept that an increase in 
the number of sensors does not lead to improvement for quadratic surfaces 
we used 49 sensor points. Here the result was still another plane surface 
with a result similar to that for the n = 25 case. Additional runs were 
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Ratio Surface for "Symmetric" 
Atmosphere with Nine Sensors 
V-10 km 


FIGURE 16b 
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made, this time using a visual range of 23 km. The results are depicted 
in Figs. 19 and 20. As before, the number of points cluster around a 
ratio of one. 

Case II: Two Point Atmosphere 

We now consider a different atmospheric state, this time characterized 

by one large peak and a large valley at points within a frame such that 

the symmetry is broken. An illustration of this case is given in Fig. 21b 

where a large peak is in the upper left section of the frame and a valley 

is in the lower right section. The corresponding radiance surface is 

shown in Fig. 21a. Again, selecting 1, 5, 9, and 25 sensor points and 

calculating a regression surface, we determined the calculated radiance- 

to-actual radiance ratio. This is plotted in Figs. 22 and 23 for a visual 

range of 10 km and in Figs. 24 and 25 for a visual range of 23 km. Here 

we see a definite improvement as the number of sensors increases because 

the regression surface more nearly matches the actual surface. In this 

case that occurs because we have less symmetry than in the four point 

atmosphere. This is also obvious from the contour plots of Figs. 26a, 

26b, 27a, 27b, and 28. For n = 1 there is a plane regression surface 

and hence the ratios follow the contour of the x surface. For n = 5, 9, 

o 

and 25 however, the surface is approaching unity as we are finding a better 
approximation to the actual surface. 

Hence, what has been demonstrated so far is that there is considerable 
improvement in the determination of intrinsic radiances if enough sensor 
points are included in a scene and if there is a lack of symmetry. Lack 
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FIGURE 20 

Histograms of Lj. (calculated) /Lj (actual) 
for a "Symmetric" Atmosphere. 
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FIGURE 23 

Histograms of (calculated) /Lj. (actual) 

for an "Asyirmetric" Atmosphere 

Visual Range “ 10 km 
Surface Albedo * 0.15 
Wavelength - 0.55 pm 











































FIGURE 26b. Ratio Surface for "Asymmetric" 
Atmosphere with Nine Sensors 
V-10 km 
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of symmetry may not be required, however, if we can use higher order 
regression surfaces. In general, we should be able to improve the 
contours by using a very high order surface but there remains the 
problem of using that surface outside the convex region bounded by the 
sensors . 

Also, we need to investigate a spatial variation of the surface 
albedo along with an atmospheric variation. The change in u o » the single- 
scattering albedo should also be studied. 

For rapid computer processing of actual data we need a more efficient 
method than the evaluation of R and S correction functions for each pixel. 
Therefore, instead of determining t q , p, and oo o surfaces and then cal- 
culating the R and S functions we can determine irradiance, path radiance, 
and total radiance surfaces by regression analysis using sensor sample 
points. This should decrease the processing time considerably. 
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3.5 CLASSIFICATION ACCURACY 

In this section we shall determine the classification accuracy for 
a realistic scene with variable atmospheres. We shall first consider the 
effect of a horizontally inhomogeneous atmosphere on the classification 
accuracy when using typical spectral signatures. Actually, there are two 
effects attributable to the atmosphere which can cause a change in the 
classification accuracy and they may or may not be separate. If the 
atmosphere is constant in the horizontal plane then there can only be a 
mixing of the signatures as a result of surface radiation from elements 
outside the sensor's field of view being scattered by the atmosphere into 
the field of view. This effect which we shall call the interference effect 
will be treated in considerable detail in the next section. 

In the more general case of a variable atmosphere the radiance may 
change throughout the entire scene as a result of the interference effect 
mentioned above and also because the atmosphere over a given target is 
quite different from point to point within the scene. For purposes of 
clarity we shall neglect the interference part of the atmosphere here and 
consider variations due to horizontal changes in the atmospheric state. 

A spatially constant atmosphere produces no change in the spectral 
transmittance, spectral irradiance, and the spectral path radiance, i.e., 

L(x,y,z,p,<f>) = L o (x,y,p,<t>) T(z,p) + L p (z,y,4>) (32) 

where L Q (x,y ,P,<J>) is the spectral radiance of the target at the surface 
point (x,y); T(z,jj) is the spectral transmittance from the surface to the 
altitude z; and L p (z,y,<f>) is the spectral path radiance at altitude z. 
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The view angles are 0(=cos , measured from the nadir, and <f>, the azimuth 
angle between the view plane and the solar plane. Thus, if we neglect the 
atmospheric interference effect then the only effect of a constant atmosphere 
is to decrease the intrinsic target radiance and increase the path radiance. 
For a spatially variable atmosphere we have for the sensor radiance 

L(x,y,z,y,4>) * L o (x,y,vi,<{>) T(x,y ,z,y ,<f>) + L p (x,y ,s,y ,4>) (33) 

where the transmittance and path radiance are now dependent upon the position 
in the plane. The surface target radiance L Q (x,y ,y ,4>) will vary not only 
because of the intrinsic change in surface materials with distance but also 
because the irradiance falling upon the target is different due to the over- 
lying variable atmosphere. For a wide range of atmospheric states, however, 
the irradiance changes very little and we will assume that it is constant 
throughout the entire scene. 

We shall n'jw consider some specific examples of the effect which variable 
atmospheres have on classification accuracy. The accuracy of classification 
can be measured by comparing actual data at a given point in a scene with a 
standard distribution obtained by generating the means, variances, and 
covariances of real data for known materials or classes. Real data were 
obtained from soybean and corn fields in Livingston Co., Illinois for an 
atmosphere which had an estimated visual range of 23 km. The means and 
standard deviations were determined and 25 data points were selected at 
random from the Gaussian distribution for each simulated field in the scene. 
The field configuration is illustrated in Fig. 29. In field Cl some of the 
sample points were closer to the mean value of soybeans than for com and 
these were classified as soybeans. Most of the points (86%), however were 
classified as corn. 
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Cl 51 C2 52 C 3 53 



85.3% C 
89.3% S 

90.7% C 
89.3% S 

88.0% C 
88.0% S 

92.0% C 
85.3% S 

85.3% C 
85.3% S 

84.0% C 
85.3% S 


PATTERN OF SIMULATED CORN AND SOYBEAN FIELDS 


Case 1: Uniform atmosphere, V « 23 km. 

87.6% of corn and 87.1% of soybeans recognized correctly. 


FIGURE 29 
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The next case simulates a variable atmosphere in which the visual range 
changes from 23 km at the top of the scene to 13 km at the bottom of the 
scene, (Fig. 30). Each data point was then altered according to Eq. 33 

and new classification accuracies were determined. This amounts to using the 
top row of fields (V=23 km) as training fields in order to determine the 
signatures. This time the percentage of classification changed slightly 
‘ Toughout the scene and in some cases actually increased. We then used 
signatures obtained from the entire scene instead of from the top row of fields. 
The classification accuracy is shown in Fig. 31 where we see little change 
from the previous case. 

We then used an atmosphere in which the visual range changed from 
23 km at the top to 8 km at the bottom. The same procedure was used as 
before and the resulting classification accuracies are shown in Figs. 32 
and 33. Here we see a slightly greater change in the percent classifi- 
cation accuracy. 

All these results can be summarized by a histogram plot of the number of 
fields with a given classification accuracy. This is illustrated in Fig. 34 
for the three cases considered. It is clear that the distribution peaks more 
strongly as we go to clearer atmospheres, i.e., the recognition decreases 
for hazy, variable atmospheres. 

It should be realized that no correlation was used in the above analysis, 
i.e., it was assumed that the soybean and corn fields were statistically 
independent. Also, a steeper atmospheric gradient could exist in reality and 
absorption by particulates should be considered. Also, as we shall see in 
the next section the interference between scenes caused by the scattering of 

radiation can be significant. This effect has not been included in the above 
analysis. Although the change in irradiance is small, the corresponding change 
in classification accuracy if this effect were included could be significant. 
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PATTERN OF SIMULATED CORN AND SOYBEAN FIELDS 


Case 2a: Non-uniform atmosphere, V * 23 km (at top) to 13 km (at bottom). 

86.2% of corn and 89.1% of soybeans recognized correctly. 
Signatures obtained from V • 23 km data set. 


FIGURE 30 
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Case 2b; Non-uniform atmosphere, V * 23 km (at top) to 13 km (at bottom). 
87.6% of corn and 87.8% of soybeans recognired correctly. 
Signatures obtained from this data set. 


FIGURE 31 
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PATTERN OF SIMULATED CORN AND SOYBEAN FIELDS 


Case 3a: Non-uniform atmosphere, V * 23 km (at top) to 8 km (at bottom) 

83.6% of corn and 90.9% of soybeans recognized correctly. 
Signatures obtained from V = 23 km data set. 
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PATTERS OF SIMULATED CORN AND SOYBEAN FIELDS 


Case 3b: Non-uniform atmosphere, V * 23 km (at top) to 8 km (at bottom), 

87.6% of corn and 87.8% of soybeans recognized correctly. 
Signatures obtained from this data set. 


FIGURE 33 
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Case 2: Non-uniform atmosphere, 13 <. V <. 23 km 



Case 3: Non-uniform atmosphere, 8 <L V <_ 23 km 

Mean recognition =■ 83.6% 



NUMBER OF. FIELDS VERSUS CLASSIFICATION ACCURACY 

FIGURE 34 
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4 

TARGET-BACKGROUND INTERACTION 

In this section we shall consider an atmospheric effect which under 
certain circumstances can be more important in multispectral data pro- 
cessing than the spatial variation of the atmosphere treated in the 

previous section. It has been noted, Nalepka et al. fe ] , that there is 
an influence of the background on the target. Such an effect is to be 
expected, particularly for hazy atmospheric conditions since radiation 
reflected from surfaces outside of the instantaneous field of view will 
eventually find its way into the sensor as a result of atmospheric scatter- 
ing. In this section we are concerned with the quantitative determination 
of this interaction effect. 

4.1 THE GENERAL EQUATION 

In a mathematical analysis of a problem as complex as the transfer 
of radiation through a medium with inhomogeneous boundary conditions it 
is important that one analyze carefully the mathematical procedures to be 
used so as to understand the limitations of a model based upon these 
calculations . 

The general equation which is used to describe the spectral radiance 
at a sensor is 

L = L I + L„ (34) 

o P 

where L is the radiance at the target in the direction of the sensor, T 
o 

is the transmittance between the target and sensor, and L p is the path 

l 8 ^R. F. Nalepka, H. M. Horwitz, N. S. Thomson, Report No. 31650-30-T, 
Willow Run Laboratories, The University of Michigan, Ann Arbor, 1971. 
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radiance, which is the result of radiation having been scattered by the 
atmosphere into the direction of the sensor. In a detailed analysis of 
path radiance as we shall undertake in this section it is desirable to 
subdivide this term into additional parts, i.e., 

L p - L p + L p + L p (35) 

o t s 

where Lp is that component which represents all multiply scattered 
o 

radiation which has never reached the surface, L p is the component 
that represents radiation originating (by reflection) from the target 

and being scattered by the atmosphere one or more times. The average of this 

component over the solid angle of the sensor can be considered as a part of 

the target radiance. Finally, Lp is the path radiance component arising from 

s 

radiation originating from the background elements and being scattered one or 
more times. These components are illustrated in Fig. 35 for the case of 
single scattering. 


Sensor 



FIGURE 35. RADIANCE COMPONENTS IN A SCATTERING ATMOSPHERE 
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The radiance at the surface, L q , is not independent of the atmosphere. 
In fact, it depends on the overlying medium in the following way 

V») = j 5 [ L S olat ( ®' ) + I 'aky^ , 'j dS ’ (36) 

^ A f 

where p(ft,ft ) is the bidirectional reflectance of the surface and 

A A A ^ 

Lsolar^ anc * are t * ie so ^ ar an< * s ky radiances at the ground, ft 

A 

and ft are the incoming and outgoing directions of the radiation. If we 
have a Lambertian (perfectly diffuse) surface then the surface radiance 
is given by 



(37) 


where p is the direction independent hemispherical reflectance (albedo) 

and E is the total do .welling spectral irradiance on the target. 

The pure atmospheric path radiance L p and the irradiance E will be 

o 

calculated according to the usual radiative-transfer model. The primary 

dependence of the total radiance L on an inhomogeneous surface arises 

- 

from path radiance terms L p and L p . It is the determination of these 

t s 

components which will be the object of our mathematical analysis. 

4.2 MULTI-DIMENSIONAL RADIATIVE-TRANSFER 

If the atmosphere has spatial variations but has a uniform surface 
we must have an inhomogeneous equation with inhomogeneous boundary con- 
ditions. On the other hand, for a uniform atmosphere and a non-uniform 
surface we also have an inhomogeneous equation with inhomogeneous boundary 
conditions. We shall specialize our analysis to consider a time-independent, 
spatially uniform medium in the horizontal plane bounded by a non-u' >Iform 
Lambertian surface. 
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4.2.1 THE TRANSFER EQUATION 


We shall consider a transport equation 

HL(r,0,t) - Q(r,n,t) (38) 

where L(r,fi,t) is the spectral radiance at position r with direction 0 at 
time t and Q(r,fi,t). is the source of radiation. H is a lir.ear Boltzmann 
operator given by 


H E c It + “* V + K(? ^ ) 


g(r,0) ' 
4ir 


P (r,n,n') dn* 


(39) 


where K(r,0) and B(r,n) are the extinction and scattering coefficients 
and p(?,fi,fi') is the phase function. We will not go into the details of 
the interesting symmetry properties of the transport equation but refer 
the reader to Case and Zwelfel [9 ] for such explanations. If we have 
time independence and an isotropic, rotationally invariant medium then 
the radiative-transfer equation is 

n-VL(r.n) + ic(r)L(f,n) - j p(r,n.n , )L(?,fj , )dfi' - Q(r,fi) (40) 

» i 
0 

where Q(r,ft) is the distributed source. For a monodirectional point 
source of unit strength we have 

Q(r,n) - 6(r - r o ) 6 (5 - fl o ) (41) 

A A f 

We consider the adjoint equation formed by replacing fi and fl with 
their time reversed states, i.e., — 

-n.VL(r,-n) + k ($)£(?, -W) - J p(?,n.ft')L(r,-fl*)dn' - Q(r,-n) (42) 

fi' 

where L(r,-n) is the adjoint cf L(r,0) . 

^K. M. Case and P. F. Zwelfel, Linear Transport Theory, Addison-Wesley 
Publishing Co., 1967. 


4.2.2 GENERAL SOLUTION 


The solution of the transport equation can be found by multiplying 
the transport equation by L(i,-ft), the adjoint equation by L(r,ft), sub- 
tracting and applying Gauss' theorem, i.e.. 


V.U dV 
V 


J 

A 


n^U dA 


( 43 ) 


where V is a volume, A is the surface area and n.^ is the unit normal vector. 

If we let 

L(r,ft) = G(r,ft; r^) (44) 


and 


Q(r,ft) = 6(r - f Q ) 6 (ft - ft Q ) 


(45) 


where G(r,ft; r o ,ft Q ) is the adjoint Green's function, then using the 
reciprocity relations for Green's functions we then find the general solution 
to be the following; 


L(r,ft) = 


G(r,ft; r ? ,ft ') Q(r’ ,ft')df ’dft' 


V ft 


G(r ,ft ; R\ft’) n ± (&*)•«' L(R* ,ft')dA dft* 


(46) 


A ft 


It should be noted that Q(r,ft) is the distributed volume source and L(R,ft) 
is the radiance on the surface. Hence, Eq. (V<>) is the general solution 

• ^ . 3 

of the radiative transfer equation of any distribution of sources within or 
on the boundary of the medium. The main problem then, is to determine the 
volume and surface Green's functions. 


71 
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-1 


4.2.3 GREEN'S FUNCTION 

The radiative-transfer equation that we must solve is the following: 


n*VL(x,y,z,y,<}>) + k(z) L(x,y,z,y,$) 


2tt 1 

~ j } •♦') L(x,y,z f y = Q(x,y,z,y,4>) 

0 -1 (47) 


where 0(=cos y) and <|> are the usual spherical coordinates of polar (or zenith) 
angle and azimuthal angle respectively. In order to find the solution 
L(x,y,z,y ,<j>) for the distributed source Q(x,y ,z ,y ,4>)^ we must first determine 
the Green’s function which is a solution of the following equation: 

f2»VG(x,y,z,y,<j>) + <(z) G(x,y,z,y ,*) 


2 it 1 

_ g_(j?l 

4tt J j 

0 -1 


p(z»y» < t>>y* » < j> , )G(x,y,z,y ' .^^dy'd^' 


6(x - x ) 6 (y - y ) 5(z) 6(y - y ) 6(<{> - <j> ) 
o o o o 


where (x ,y , 0) are the coordinates of the surface source and 0 (=cos y ) 
o y o o o 

and <Jj^ are the corresponding angles for the monodirectional source. Since 
we have homogeneity in the x-y plane we can use the two-dimensional Fourier 


transform, i.e.. 


G(k,£,z,y ) = | | e i ^ kx+ ^ y ^G(x,y,z,y,4>)dxdy 


and an inverse. 


G(x,y,z,y,<f>) = — f [ e i ^ lcx+ ^ y ^G(k,£,z,y ,«}))dkd£ (50) 

(2") 1 1 
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Taking the Fourier transform of Eq. (*(5?) and finding the formal solution 
we get for the general Green’s function transform 


G(k t Z , z ,u » 4*) 


ri rri 

1-y ^ r -- ^ M zsin<t>]£ 


e -[t -t(z)]/m i[ x Q + zcos^jk i[y o + 

— e e 

y 


6(y-u Q ) 6($-<t> o ) 


z2f 1 


4iry 



» 


J 




0 0-1 


dp d<j> 'dz ' 


(51) 


where /i-y^[k cos<j> + & sin<J>], 

4.2.4 SINGLE-SCATTERING 

Equation (<$•) is an integral equation for the transformed Green's 
function. It represents a complete solution with an infinite number of 
scatterings for a monodirectional point source. Iterating once gives us 
the single-scattering solution for the transformed Green's function. Taking 
the inverse transform then gives us the solution, i,e.. 


G s (x,y,z, U ,*i VVVV = 


1 

4-.fy 


z 

, p(z* ,y,ij>,y 0 ,4> 0 ) -[T(z’)-T(z)]/y -[T o -T(z')]/y o 

8(z ) — e e 

M o 

0 


x 6 


/l-M 2 

o * . . /l-y , 

x-x z cos<}> (z-z )cos<J> 

o y o y 

_ o 


J ^ 

/l-y , », 


£7 


o » 

y-y z cos4> 

J J o y o 


,, 1 V, ' U (z-z')sin<{> 


i 

dz * 


(52) 
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The geometry appropriate to Eq. (5*-) is Illustrated in Fig. 36. An 

observer o'* sensor is located at point (x,y,z) and receiving direct radiation 

from the point (x ,y ,0) on the surface. h«*diation also emanates from either 
s s 

a target or background point source designated by (x Q ,y o ,0), is attenuated 
to point (x’.y’.z’), and is then scattered into the direction (6,<j>) and is 
attenuated to point (x,y,z) , 

4.2.5 SURFACE PATH RADIANCE 

Since we have no volume source of radiation we consider our only source 
to be the surface. Hence, the solution is 

L p (x,y,z,vi,<J>) = 
s 


2 IT 1 

I J J I wvVo L <vvv*o> 4 vW y o (53) 

*o y o 0 0 


Inserting the singly-scattering Green's function and carrying out the inte- 
gration over the angular coordinates we find the path radiance arising from 
the surface to be 

z 

_ . . 1 *3 B(z ) , » % 

L p (x,y,z,u,<(>) - £inT y t2 p ^ z » cosx ^ 

s z 

• 

x y 0 
o ■'o 

-[t(z')-t(z)]/m -[t 0 -t(z')]/m* * * , 

x e e L(x ,y ,4> )dz dx dy^ (54) 

o o o o 
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FIGURE 36. GEOMETRY FOR SINGLY SCATTERED SURFACE RADIATION 
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The angle x is the scattering angle. 

Equation (54) is therefore a mathematically exact formula for that 
portion of the path radiance which arises from radiation reflected from the 
surface and is single-scattered . The total spectral radiance at a sensor is 
then given by 


L(x,y»z,w,t) = L q (x ,y g ,u,<fi)T(z,w) + L p (z,u,<|>) + L p (x,y ,z,u ,♦) (57) 

o s 

where L_ (z,M,<f>) is the pure atmospheric path radiance and L (x ,y ,y,4>) is 

£ OSS 

o 

the surface radiance. 

Before we analyze the effect of background on target via atmospheric 

scattering we shall look at another effect associated with the finite 

aperture of a sensor. In order to determine the mean radiance received by 

a sensor with a finite opening angle 0 we must integrate the total radiance 

s 

over that angle, weighted by the cosine, i.e., 


2 it 1 


0 

2 IT 1 


M L(vj,<J>)dMd<J> 


u dpd<j> 


J J 

0 u 


( 58 ) 


$ 

« 
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The question arises, if 6 (=cos ) is measured in milliradians how much 

s s 

does L(y,4>) differ from L ? If part of the instantaneous field of view 
lies near a highly reflecting area then it is possible that a significant 
difference does occur. The various components of interest can be illustrated 
as in Figure 37. 



FIGURE 37. RADIANCE COMPONENTS FOR A FINITE APERTURE 
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In a previous study. Turner [2], it was assumed that all radiation from 
the surface was scattered along the direct ray from the center of a target 
to the sensor. This corresponds to components 1, 2, and 3 in Fig. (37). 

The first component is the directly attenuated radiation from the surface, 
the second is the singly-scattered radiation from the target, and the third 
is the singly-scattered radiation from the background. Likewise components 
4, 5, and 6 are the corresponding components for radiation being singly- 
scattered within a erne. For hazy atmospheres the radiation is strongly 
peaked in the forward direction and therefore the radiance decreases quite 
rapidly as a function of angle. If we evaluate the radial integrand of the 
path radiance we find that it decreases as in Fig. 38. It should be noted 
that the largest part of the integrand lies within 80 cm or so of the 
instantaneous point of view. Hence, for a resolution element of ^80 m on a 
side, corresponding to the LANDSAT field of view we see that the edge effect 
is very small. 

Although the value of the integrand decreases rapidly with increasing 
angle, the value of the total integral or the path radiance will increase 
due to the fact that more surface area is being included as the distance 
increases from the center. This effect is illustrated in Fig. (39) for three 
different atmospheres. It should be noted that 60% of the contribution is 
reached at a distance of 0.8 km or ten LANDSAT pixels for a hazy atmosphere 
whereas 60% of the contribution occurs at four LANLSAT pixels for a very hazy 
atmosphere. Thus, although the total effect is less for a clear atmosphere 
than for a hazy atmosphere it is necessary to go out to greater distances to 
obtain the same percentage contribution. This is due to the fact that the 

^R. E. Turner, Radiative Transfer in Real Atmospheres, Report No. 
190100-24-T, Environmental Research Institute of Michigan, Ann Arbor, 1974. 
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RADIAL INTEGRAND 












FIGURE 39 

FRACTION OF THE SURFACE PATH RADIANCE 
AS A FUNCTION OF DISTANCE FROM ORIGIN 
FOR A UNIFORM SURFACE REFLECTANCE. 
VISUAL RANGES ARE 3 KM, 10 KM AND AO KM 
WAVELENGTH =0.55 pm. 








probability for scattering is less for a clear atmosphere but the probability 
that the radiation is scattered at a large angle is greater because of the 
smaller degree of anisotropy in the single-scattering phase function. This 
effect is illustrated in Fig. 40. Hence, the total effect of a background 


Isotropic and Rayleigh Functions 
are Multiplied by Ten 


30° 60° 90° 120° 

330° 300° 270° 240° 



FIGURE 40 ANGULAR DEPENDENCE OF SINGLE -SCATTERING PHASE FUNCTIONS IN 

ANY AZIMUTHAL PLANE 


is important and we should now consider a more realistic surface pattern 
characteristic of agricultural areas . 

4.3 SIMULATION OF SATELLITE DATA 

We shall now consider surface conditions which are more realistic than 
those used in Section 3. There we assumed a uniform surface with a background 
all composed of the same material. In this section we apply the mathematical 
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technique^ of Green's functions in order to simulate satellite multispectral 
data appropriate to extended agricultural areas . 

4.3.1 REFLECTANCE PATTERNS 

Using the mathematical methods of the preceding section we can calculate 
the singly-scattered radiance from any complex pattern of reflectance elements. 
For computational simplicity, however, we shall consider a simple checkerboard 
pattern of square fields with the target being the center square. The radiance 
in a vertical direction at space altitudes is given by 

L(x f ,y f ,<j>) = L Q (x f ,y f ><£)T 

+ L p (x f ,y f ,<|>) + L p (x f ,y f ,<J>) (59) 

o s 

where (x r t y f ) is a point within the target square , For homogeneous atmospheres 

L P = L p 

o o 

The surface path radiance, L p is 

s 


L P 

s 



( x f »y £ ! 


•V y o> 


dx dy 
o J o 


(61) 


where we must integrate over the entire surface, ( x 0 »y Q ) being any point 
on the surface. The quantity F(x^,y£,4 ! »x o ,y o ) is given by 


F(x f ,y f *.x o ,y 0 ) 


1 _ 

4tt 


,'2 


-[T( Z ')-T(Z)] -[T rt -T(z')]/y* 

e e 




(62) 
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J 


■HMumm*-. 


L 


where 


U" = — 


( 63 ) 


Kx f -x o ) 2 + <y f -y 0 ) 2 + '*' 2 


* -l 
4> = sin 



y f - y 0 


( 64 ) 


) 2 + (y f -y Q ) 2 + z ' 2 


In the case of Lambertian surfaces 


P ii y s E 


( 65 ) 


where p,. is the reflectance of the (i,j) square, u is the cosine of the solar 
ij s 

zenith angle, and E is the irradiance. The surface path radiance is therefore 
given by 


N N 

h <♦» - I I Vi, 

s i=l j=l J J 


( 66 ) 


when 1^ is the integral over the (i,j) square. It should be noted that 
the integrals I „ are independent of reflectance so that the lengthy compu- 
tations for these integrals need only be done once for a specific set of 
wavelengths, atmospheric states, and field sizes. Any number of reflectance 
patterns may then be considered by multiplying each by each 1^ . If we 
look at the origin of coordinates, i.e., at the center of the target square, 
then by symmetry 


h s ■ p oo I oo + I’lJ + p j ,i + P l,-j + P -J ,i + p -i,-j + P -J ,-i < 67 > 


” (Pj 4 + P^ 4 + P'j 4 + P 4 •) (6.4 +6. +6. -6. 6 )] I 

i»j i,~3 -i.j -i»~j ij oj oi oi oj' J 


ij 
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! 1 ! ] — U. . __i 

where the 6.. are Kronecker deltas and p Is the reflectance of the target, 
ij oo 

It should be noted that only (N-l)(N+5)/8 integrations need be performed 
for an N by N grid. 

We considered three atmospheric states characterized by visual ranges 
of 3 km, 10 km, and 40 km corresponding to very hazy, hazy, and relatively 
clear conditions respectively. The wavelengths chosen are those corresponding 

« 

to the centers of the four ERTS channels, i.e., 0.55 ym, 0.65 ym, 0.75 ym, 
and 0.95 ym. It is assumed that the sensor is looking in the nadir direction 

from space at fields in the midwestern part of the United States on 24 June 
1975. The atmosphere is uniform over the extended area and contains no 
absorbing particulates. We chose, for computational convenience, a 19 by 19 
grid, i.e., 361 square fields, each 2 pixels or 160 m on a side. Thus, each 
field contains 4 resolution elements. Configuration 1 i3 depicted in Fig. 41. 

We consider a target in the center of the grid and a variable background. 

A new reflectance is imagined as moving in from the right, column-by-column 
unti. . completely overtakes the original background. Configuration 2 is 
depicted in Fig, 42. Here we consider a new background which is expanding 
outward from the target field, ring-by-ring until the original background is 
completely overtaken. Any configuration is possible but we shall consider 
these two examples. 

First, we analyze a black target and a variable black-white background. 

Figure 4.3 illustrates the variation in the ratio of the surface path radiance 
to the pure atmospheric path radiance as a white cover advances over the 
scene. It is assumed, however, that the target Is always black. It should 
be noted that there is a very large increase only after the white cover 
passes the target and then approaches an asymutctic value. Figure 44 illus- 
trates the spectral variation of the same ratio. It is almost linear with 
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FIGURE 43 

RATIO OF SURFACE PATH RADIANCE TO 
ATMOSPHERIC PATH RADIANCE AS A 
FUNCTION OF A VARIABLE BLACK-WHITE 
BACKGROUND AND A BLACK TARGET; 
WAVELENGTH =0.55 pm, SOLAR ZENITH 
ANGLE * 31°, VISUAL RANGE - 10 KM. 
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FIGURE 44 

RATIO OF SURFACE PATH RADIANCE TO 
ATMOSPHERIC PATH RADIANCE AS A 
FUNCTION OF WAVELENGTH FOR A 
VARIABLE BLACK-WHITE BACKGROUND 
AND A BLACK TARGET; SOLAR ZENITH 
ANGLE = 31°, VISUAL RANGE = 10 KM 
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wavelength but the slope of the curves is quite different depending on the 
amount of white background. 

Second, we shall consider the same configuration only this time we have 
a variable wheat-grass background where reflectances for healthy, mature wheat and 
grass were used. This would simulate a scan across a wheat-grass boundary. 

Figure 45 illustrates the effect. Pure wheat means that the entire scene is 
wheat. Edge of wheat means that we are looking at a wheat field target with 
one half of the scene covered by wheat and the other half by grass . As we 
cross over the boundary to the grass one can see that a large change occurs 
due to the dominance of the radiance from the target itself. Then as we 
progress further into the grass region the radiance changes slightly. 

Third, we consider a wheat target being surrounded by rings of grass. 

Figure 46 depicts the variation in the spectral radiance in this situation. 

We have a wheat field surrounded by a black background. Then we have an 
all wheat scene which clearly illustrates the magnitude of the background 
interference. It should be noted that this represents about 90% of the 
radiance possible, i.e., an infinite plane of wheat would raise the values 
by 10%. The third curve represents a wheat field surrounded by a total 
grass scene. The percentage change is significant but this will have to 
be translated into actual signatures before classification accuracies can 
be determined. 

Fourth, it is perhaps more useful to consider the per cent change in 
radiance from one configuration to another. Let us imagine a wheat target 
surrounded by a grass background. Figure 47 shows that there is about a 
25% change in channel 3. Now, if we have a wheat target and the next ring 
of elements (fields) is also wheat with the rest of the surface being grass, 
then the change is only M.7% in channel 3. It continues to decrease, of 
course, as we add more wheat to the area surrounding the target. 
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FIGURE 45 

TOTAL SPECTRAL RADIANCE AS A FUNCTION 
OF A VARIABLE WHEAT-GRASS BACKGROUND 
FOR A WHEAT TARGET AND A GRASS TARGET 
VISUAL RANGE - 10 KM, SOLAR ZENITH 
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FIGURE 46 

VARIATION OF WHEAT RADIANCE FOR A 
BLACK, WHEAT, AND GRASS BACKGROUND 
VISUAL RANGE - 10 KM, SOLAR ZENITH 
ANGLE = 31°. 
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FIGURE 47 

PER CENT CHANGE IN TOTAL RADLANC 
FOR A WHEAT TARGET SURROUNDED BY 
RINGS OF WHEAT AND GRASS. 

SOLAR ZENITH ANGLE = 30°, 

VISUAL RANGE - 10 km 
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How much does the per cent change vary as a function of sun angle? 

This is illustrated in Fig, 48 for four solar zenith angles. There seems 
to be a maximum effect in the range 6 o ^ 30°-60°, which is also the range 

r 

of sun angles for which most data are collected. 

In order to see the simultaneous change in total radiance and surface 
path radiance we can plot one versus the other for all four channels. This 
is depicced in Fig. 49. We go from an all wheat area to a wheat field 
surrounded by grass. The length of the displacement depends upon visual range 
the greater the visual range, the smaller the line segment. Thus, for no 
atmosphere the line segment is zero . 

The greater the reflectance of the background relative to the target, 
the greater is the change in total radiance. Thus, for a dark target 
such as soil or loam the change can be quite large . This is shown in 
Fig. 50 for soil being surrounded by various amounts of grass . A corre- 
sponding vs. Lp plot is depicted in Fig. 51. 

s 

In order to see the per cent change for various materials we can analyze 
wheat targets with grass backgrounds, soil targets with grass backgrounds, 
and diseased wheat targets with a wheat background. This dramatic effect 
is illustrated in Fig. 52. 

What is the effect of different atmospheric states on the radiance 
values? In Fig. 53 we consider soil targets with various backgrounds. As 
is to be expected, the change increases with increasing turbidity, or 
decreasing visual range. 
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FIGURE 48 

PER CENT CHANGE IN RADIANCE 
BETWEEN A VfHEAT FIELD 
SURROUNDED BY GRASS AND 
AN ALL WHEiT AREA FOR FOUR 
SOLAR ZENITH ANGLES. 

VISUAL RANGE = 10 KM 


WAVELENGTH (ym) 
















RADIANCE (mW/cm -ster-ym) 



P!S0lHfii5»3BI8Ha8E! jijft B . ij jtij 




SOLAR ZENITH ANGLE = 30 


SURFACE PATH RADIANCE (mW/cm -ster-ym) 


FIGURE 49 


TOTAL RADIANCE VS. SURFACE PATH RADIANCE 
FOR WHEAT-WHEAT TO WHEAT-GRASS TRANSITION 
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(shortest vector) 
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FIGURE 50 

CHANGE IN LOAM RADIANCE CAUSED BY 
SURROUNDING GRASS. TARGET AND 
RING WIDTH = 160 METERS, VISUAL 
RANGE = 10 KM, SOLAR ZENITH 
ANCLE = 31°. 
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FIGURE 51 

L_ VS . L_ FOR A SOIL TARGET 


SURROUNDED BY RINGS OF GRASS 
SOLAR ZENITH ANGLE = 30°, 
VISUAL RANGE = 10 km 


ALONE 


L_ (mW/ cm -ster-um) 
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FIGURE 52 

PER CENT CHANGE IN RADIANCE FOR 
SOIL SURROUNDED BY GRASS, 
DISEASED WHEAT SURROUNDED BY 
WHEAT, AND WHEAT SURROUNDED BY 
CRASS. SOLAR ZENITH ANGLE = 30' 
VISUAL RANGE = 10 KM 
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FIGURE 53 

VARIATION OF LOAM TARGET WITH 
VISIBILITY FOR A BLACK, LOAM, 
AND GRASS BACKGROUND. WAVE- 
LENGTH =0.75 ym, SOLAR ZENITH 
ANGLE = 31°. 
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Fin., ^y, we should consider the various components which make up the 

total racHancu. Figure 54 illustrates the spectral variation of the beam 

radiance L , the pu. * at^o-oheric path radiance L p , the surface path 
B o 

radiance L p , and the tOw t radiance L^. Although the dominant term is 

the beam radiance, the surface path radiance is by no means insignificant. 

For shorter wavelengths or smaller visual ranges the multiply scattered 

surface path radiance can be expected to increase. Thus, the background can 

certainly alter the total radiance. The pure atmospheric path radiance is 

independent oi surface conditions and decreases monotonically with increasing 

wavelength . 
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FIGURE 54 

RADIANCE COMPONENTS FOR A WHEAT 
TARGET, GRASS BACKGROUND. 

SOLAR ZENITH ANGLE = 30°, 

VISUAL RANGE - 10 km 















5 

INTERACTION CORRECTION ALGORITHM 

In this section we shall develop mathematical techniques which can 
be used to correct actual muitispectral data for the interaction effect. 

It is assumed that the atmosphere is homogeneous in the horizontal plane 
but the flat surface may have any degree of spatial inhomogeneity. 

5.1 GENERAL SURFACES 

We shall first consider general, non-Lambertian surfaces. The general 
equation is 

L(r,fl) = L o (r s ,fl)T(r,S2) + L p (r, 0 ) + L p (r,Q) (68) 

o s 

For homogeneous atmospheres 


T(r,fl) = T(z,fi) 

(69) 

L p (r ,0) = L p (z,0) 

(70) 


o o 


Thus , 

L(r ,z,0) = L (r s ,6)T(z,fi) + L p (z,fi) + L p (r ,z,6) (71) 

p o s F 

where the vectors r and r are illustrated in Fig. (55). Now, the singly- 

p s 

scattered path radiance is given by 


h p ( r »z,f2) 
s v 


J J 


G(r ,z,0; r ,0 ) L (r ,0 )dr aO 
o p* * o o o o o o o 


(72) 


Taking the Fourier transform of Eq. (71) leads to the following equation: 

-ik 


L(k,z,0) * T(z,0) e 


t * 

L o (k,fl) 


+ (2ir) z 6(k) L p (z.fl) + L p (k,z,0) 


(73) 
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. GEOMETRY FOR CORRECTION ALGORITHM. 

Sensor Location at r. 

Ground Observation Point at r . 

s 


FIGURE 55 


where the bar indicates the Fourier transform. Now, because of the trans- 
lational invariance of the Green's function we get 

L p (k,z,0) - 
s 


■ 

n*fi G(k,z,fl,0 ) L (k,n )dft (74) 

o o o o o 


Hence, we finally get an integral equation for the surface radiance transform, 
i »e • , 


r ,, s , itt-r <2 '> S(k) h- (2 ' a)e 

l (k,n) - e t ...... 

T(z,sj) I(z M 


i£*r. 


T(z,fl) 


n*fl G(k,z,R,ft ) L (k,S2 )dfl 
o o o o o 


(75) 


Taking the inverse Fourier transform leads to the following result for the 
corrected radiance: 


L (r ,0) 
o p* 


L(r u ,z,fi) - L p (z,0) 
o 

T(z,n) 


n*n G(r ,z,fi,0 ) L (r - r)dS2 dr' (76) 

o u * o o u u o u 

where r^ * r^ - r^ . Equation (76) says that we must subtract the pure 
atmospheric path radiance from the actual data value and divide by the trans- 
mittance and then subtract a modified surface radiance from integrated over 
the entire surface. We have an integral equation for the corrected radiance 

a 

L(r,0) which can be solved in various ways. In particular, if we use the 
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single-scattering Green ss function in Eq. (76) then a definite (but complicated) 
analytical result can be found for the corrected radiance. 

5.2 LAMBERTIAN SURFACES 


If we assume a perfectly diffuse (Lambertian) surface considerable 
simplification occurs- in the mathematical method, 'Equation (76) then 
becomes . 


A 

L o< fc > 


t’Vt' o\ ik *r 
L(k,z ,£i) t 

T(z,fl) 


o « „ 

(2tt) Z 6(1*.) Lp (z,8) e 

o 

T (z,f5) 


ik*r _ __ 

e g(k,s) h o (fc) 

T(z,8) 


(77) 


where g(k,z) is the transform of an isotropic Green’s function integrated 
over a hemisphere. Equation (77) can therefore he solved for L Q (k), i.e.„ 


L o 0O = 


L(k,z,£?) - (2ir) 2 o(k) L 


o 


ik*r. 


i£*r. 


T(z,8) + g(k 9 z) e 


(7S) 


Thus, knowing the surface pattern and the optical thickness or visibility 
conditions one can find the transmittance, pure atmospheric path radiance, 
and the Green’s function. Inverting the above transform then gives us the 
surface radiance free of the interaction effects. The general algorithm Is 
illustrated in Fig. (56). 
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CORRECTED RADIANCE 
(FOURIER TRANSFORM) 


CORRECTED RADIANCE 
(INVERSE TRANSFORM) 


FIGURE 56. ALGORITHM FOR CORRECTION OF MDLTISPECTRAL DATA 
FOR ATMOSPHERIC INTERACTION EFFECTS 
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CONCLUSIONS AND RECOMMENDATIONS 


In the current analysis wa have analyzed two atmospheric effects, 
the spatial variations in the atmosphere and the target-background inter- 
action, Of all possible atmospheric variations these will cause the 
greatest change in classification accuracy. 

The relationship between the various radiometric quantities which 
we considered in some detail is useful for the determination of atmospheric 
and environmental parameters which, in turn, are used in an atmospheric 
correction algorithm. Thus, a measurement of irradiance or sky radiance 
along vith transmittance will allow certain parameters to be estimated. 

For variable atmospheres it is necessary to perform several measurements 
at selected points in a scene. As we have seen in the simulation analysis, 
the greater the number of auxiliary measurement sites the better is the 
agreement between actual and calculated radiances. 

In the classification problem we showed that rather severe atmospheric 
optical thickness gradients are needed in order to realize a significant 
change in the per cent classification accuracy. It should be realized 
however, that no correlation was assumed between classes and hence the 
soybean and corn fields were statistically independent. Also, no absorption 
was included in the simulations made. If remote sensing data are collected 
over a large area or near an urban complex then there could very well be 
large spatial changes in the absorptive part of the optical thickness. At 
first, it seems that strong gradients in the optical thickness surface are 
unrealistic, but that is not the case if absorption is included. Also, it 
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was assumed throughout that there was no interference arising between 
fields as a result of atmospheric scattering. As was learned in the 
next section such a neglect can affect classification accuracy. Hence, 
we conclude that if multispectral data are collected over small, homogeneous 
regions away from sources of pollution and in very clear atmospheres then 
the classification accuracy is not greatly affected by atmospheric changes. 
It should he realized, however, that these are severe restrictions. A not 
inconsiderable amount of remote sensing data collection is performed under 
adverse weather conditions. Also, it is becoming increasingly more diffi- 
cult to find extended areas which are completely free of aerosol absorption. 
Thus, in order to extend the capability of remote sensing to a wider range 
of atmospheric conditions it seems appropriate that the atmospheric correc- 
tion procedures outlined in Section be implemented on actual data. 

Only in recent years has the effect of surface albedo been realized 
to be of importance in multispectral data analysis. It has generally been 
assumed that spectral signatures are essentially determined by the target 
material. Background albedo was recognized as being of some importance in 
the estimation of the magnitude of path radiance but no real quantitative 
analysis was done, especially for various surface patterns of fields. As 
a result of the work done in the current study we are now able to analyze 
the effect of varying target-background geometry on the spectral radiance 
characteristics of s target material. We can consider fields of any size 
and simulate any complex pattern. We have found that the background changes 
the irradiance on the target by a very small amount, much less than the 
effect of background on path radiance. The current analysis makes use of 
few approximations. The integration is done numerically with realistic 


anisotropic phase functions and for vertical changes in the atmosphere. 

This means that, as far as the surface path radiance is concerned, we 
have an exact formulation for a realistic, vertically inhomogeneous 
atmosphere . 

The results of the calculations clearly indicate a sizable change 
in the radiance from targets of standard materials such as wheat, soil, 
and grass when various background materials are present. It should also 
be pointed out that the effects presented here are in a sense minimal 
since we cannot integrate over an infinite surface. They do represent, 
however, most of the effect due to single scattering, i.e., the effects 
are perhaps fifty to ninety per cent of the full value. Most of the 
simulations considered were for a 10 km visual range but it is apparent 
from the graphs that even, for high visibilities the effect can be important. 
Hence, we conclude that if remote sensing is performed without taking into 
consideration the geometry and type of background elements then a false 
determination of the signature may result. This is especially true for 
highly variable contiguous fields. 

For the variable atmospheric problem we recommend that additional 
variable surfaces be considered, i.e., those which portray a larger amount 
of absorption. Also, non-linear optical thickness gradients in the hori- 
zontal plane should be used which are more indicative of variable atmospheric 
conditions near urban areas. In addition, the atmospheric correction algorithm 
outlined in Section 2 should be used on real data with known atmospheric 
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Although much has been learned in the target-background interaction 
problem there are a number of things which should be done to improve the 
accuracy. First , we used mean values of reflectance for the various fields . 
Most of these data > however are for the leaves of plants rather than for the 
crops themselves. Thus, the current model should be interfaced with a 
reflectance model for the materials used. The mathematical procedures 
developed herein are versatile enough to include the bidirectional reflectance 
of any material. The interesting goniometric properties associated with 
some crops could then induce changes in the present results, which could 
be either larger or smaller than are indicated by assuming Lambertian 
surfaces . 

It would be interesting to include the interaction effect with the 
variable atmosphere and to see the total change in classification accuracy. 

In any case, for a uniform or non-uniform atmosphere the classification 
analysis procedures should be carried out for a complex pattern of inter- 
acting fields. This simulation can be done with the present analytical 
techniques . 

Finally, the itiverse problem should be analyzed, i.e. s the implementation 
of the interaction correction algorithm outlined in Section 5. 

Including all of these effects, i.e., performing correction analysis 
followed by classification procedures should result in a significantly 
improved classification accuracy of raultispectral data. 
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